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Description 

The present invention relates to a method for the detection and identification of particles in a suspension, compris- 
ing the following steps: 

5 

a. generation of acoustic signals having the form of a beam using an acoustic source; 

b directing the acoustic signals at at least one measurement volume within the suspension, the boundaries of the 
measurement volume in the axial direction with respect to the acoustic source being defined with the aid of time 
windows; 

io c reception of acoustic reflection signals produced by reflection of the acoustic signals by the particles in the at 
least one measurement volume; 

d. conversion of the acoustic reflection signals into electrical reflection signals; 

e. counting numbers of electrical reflection signals which have an amplitude in excess of a predetermined value 
and conversion thereof into numbers of particles which are larger than a certain size. 

15 

A method of this type is disclosed in British Patent 1 ,012,010, which describes a method and equipment for count- 
ing and measuring such particles, wherein acoustic samples are taken in various measurement volumes along the 
acoustic axis of the acoustic transducer in the suspension. By using suitable time windows when receiving reflected 
acoustic signals, the particles in, for example, four predetermined measurement volumes, which are each located a pre- 

20 determined distance away from the transducer, are counted. By making use of a threshold voltage which the electrical 
signals produced from the acoustic signals must exceed in order to be counted, which threshold voltage is different for 
each zone, a minimum size for the particles to be counted is selected for each zone. Assuming that the particle distri- 
bution is the same in each zone, a rough estimate of the number of particles, subdivided according to particle size, can 
be obtained using this known method and using simple mathematical methods. 

25 US Patent 3 774 717 describes a method and equipment for the detection and identification of small particles; for 

example biological cells, which, for example, are located in a medium which flows transversely to the direction of prop- 
agation of an acoustic signal. Each of the particles gives a specific scatter of the acoustic signal, depending on the size, 
the shape and the acoustic impedance of the particles. The acoustic signal has a wavelength and an effective cross- 
section of the order of magnitude of the particles to be detected and identified. Blood cells, for example, are detected 

30 with the aid of an acoustic signal of 860 MHz. Therefore, particles can be identified with the aid of techniques which are 
known from radar technology. The technique disclosed in this patent is unsuitable for in vivo detection and identification 
because the wavelengths used allow only very restricted depths of penetration in biological tissue. 

The methods described above are based on the ultrasonic pulse-echo technique. With this technique use can be 
made of a so-called ultrasonic transducer, which converts an applied electrical pulse into an acoustic (ultrasonic) signal 

35 and which is also capable of converting an acoustic (ultrasonic) signal which is incident on the reception surface back 
into an electrical signal. Therefore, the transducer serves as transmitter and as receiver for ultrasonic signals. It is also 
possible to use independent transmitters and receivers. A high-frequency pulse-echo recording of the flowing suspen- 
sion is made. A focusing transducer can be used for this. In Figure 1 the "illuminating" sound signal is shown diagram- 
matically in an arrangement known per se. 
•40 As is shown in Figure 1 . the principal axis z of the sound beam 20 generated by a transducer 23 is perpendicular 
to the direction of flow P of the suspension 21 flowing in a channel 24 and, consequently, to the direction of movement 
of the particles 22 present in said suspension. When a particle 22 passes through the sound beam 20. the incident 
sound field will be reflected by said particle and the reflected signal will be captured by the transducer 23. The received 
signal is converted by the transducer 23 into an electrical signal, which is transmitted to the transmission and reception 

45 electronics 25. The transmission and reception electronics 25 transmit the signal to a computer 26, which is connected 
to a memory 27 for storing measurement data. The computer 26 is provided with suitable software for evaluation of the 
measurement data. The electrical signal from a single measurement is indicated diagrammatically by s and is a time 
signal, the time axis t indicating the propagation time of the sound. The response of the particle 22 in the measurement 
volume can be detected in the recording at that moment in time which corresponds to the propagation time of the ultra- 

50 sonic pulse between transducer 23 and reflecting particle 22 and back again. 

Only reflections within an applied time window [t 1t t 2 ] (see Figure 1) are processed in the analysis. The measure- 
ment volume is thus limited in the axial direction by z , =\,xJ2 and z 2 = t 2 .c/2 (c is the speed of propagation of the 
sound). In the lateral direction the measurement volume is limited by the shape of the acoustic beam. 

A method for counting the number of particles using a measurement set-up of this type is described in the above- 

55 mentioned patents and in Groetsch, J.G.: "Theory and application of acoustic particle monitoring systems", Jr. 
Advances in Instrumentation and Control 45 (1990), Part 1.. Generally speaking, with this method a specific threshold 
value is chosen for the amplitude of the reflected signal. If Ihe recorded signal within the time window [t 1 , t 2 ] is in excess 
of this threshold value, this is interpreted as the pr sence of a particle in the measurement volume: On condition that 
the particle concentration is so low that the risk of the simultaneous presence of more than one particle within the meas- 
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urement volume is negligible, the particle concentration can be estimated by counting the number of recordings in 
excess of the set threshold value. 

This method takes no account of variations in particle size and no distinction is made between different types of 
particles which are possibly present in the suspension. 
5 The aim of the present invention is to provide a reliable and more precise method for the characterisation of a sus- 

pension and the particles present in the suspension. In this context characterisation is understood to be estimation of 
the particle concentration, the particle size distribution, the shape of the particles and the reflectivity of the particles. 

A further aim of the invention is to provide a method for the characterisation of suspensions and the particles con- 
tained therein, with which method account does not necessarily have to be taken of the condition that no more than one 
w particle may be present in the measurement volume at any one time. 

The first-mentioned objective is achieved with a method of the above-mentioned type, which method is character- 
ized by the following steps: 

f. composing at least one curve on the basis of a cumulative count of the number of reflection signals which have 
is an amplitude in excess of a specific value as a function of the amplitude; 

g. comparison of the at least one curve with predetermined standard cumulative count curves and deduction of at 
least one feature from a set of features comprising: material properties, particle concentration, particle shapes, par- 
ticle size and standard deviation thereof and particle size distribution. 

20 Material properties thus deduced may, e.g., be density and compressibility. 

A method of this type can be used successfully where the suspension flows relative to the acoustic source, where 
the acoustic source is moved relative to the suspension (for example along the suspension), where, for whatever rea- 
son, the properties of the suspension change as a function of time and the acoustic source is fixed, and where an array 
of acoustic sources is used instead of a single acoustic source, the acoustic sources always being activated in succes- 

25 sion. 

For the purposes of the further objective, one embodiment of the method according to the invention is character- 
ized in that: 

in step a, the acoustic beam generated has such a large aperture, and the time which elapses between two suc- 
30 cessive measurements is so short, that each particle is exposed several times while passing through the beam and 
that, depending on the lateral position of a particle in the measurement volume, a varying angle-dependent reflec- 
tion of the acoustic signal is produced; 

prior to step f, one or more different types of particles present in the suspension are identified on the basis of a 
series of angle-dependent reflection signals received successively over time; 
35 - when composing the curve in step f, the maximum value of the amplitudes of a series of successive angle<Jepend- 
ent reflection signals, received over time, from a detected particle from the one or more groups is taken as the 
amplitude of the electrical signals, which are produced after conversion of the acoustic reflection signals from that 
particle. 

40 With a method of this type, the measurement volume chosen is so large that, within the measurement volume, the 
angle of incidence varies as a function of the lateral position. If a particle in the flowing suspension is "exposed" various 
times in succession by an acoustic signal in pulse form, the successive reflection signals differ as a consequence of 
angle-dependent reflection. Said angle-dependent behaviour is highly dependent on the shape of the particle, as can 
be seen from, inter alia: M.G.M. de Kroon: "Acoustic backscatter in arteries - Measurements and modelling of arterial 

45 wall and blood", thesis 1993, ISBN 90-9006182, Section II. Therefore, a particle can be characterized by this means. 

In an alternative embodiment the present invention relates to a method for the detection and identification of parti- 
cles in a suspension, comprising the following steps: 

a. generation of acoustic signals using an acoustic source; 
so b. directing the acoustic signals at at least one measurement volume within the suspension, the boundaries of the 
measurement volume in the axial direction with respect to the acoustic source being defined with the aid of time 
windows; 

c. reception of acoustic reflection signals produced by reflection of the acoustic signals by the particles in the at 
least one measurement volume; 
55 d. conversion of the acoustic reflection signals into electrical reflection signals; 

e counting numbers of electrical reflection signals which have an amplitude in excess of a predetermined value 
and conversion thereof into numbers of particles which are larger than a certain size- 
characterized in that the method also comprises the step of applying an inversion algorithm on the ampli- 
tudes of the electrical reflection signals to deduce at least one feature from a set of features comprising: material 
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properties, particle concentration, particle shapes, particle size and standard deviation thereof and particle size 
distribution. 

The present invention also relates to equipment for the detection and identification of particles in a suspension, 
5 comprising: 

a. an acoustic source for the generation of acoustic signals in pulse form; 

b. means for directing the acoustic signals at at least one measurement volume within the suspension, the bound- 
aries of the measurement volume in the axial direction with respect to the acoustic source being defined with the 

io aid of time windows; 

c. means for receiving acoustic reflection signals produced by reflection of the acoustic signals by the particles in 
the at least one measurement volume; 

d. means for converting the acoustic reflection signals into electrical reflection signals; 

e. means for counting numbers of electrical reflection signals which have an amplitude in excess of a predeter- 
15 mined value and for converting said count into numbers of particles which are larger than a certain size; 

characterized in that the equipment also comprises: 

f. means for composing at least one curve on the basis of a cumulative count of the number of reflection signals 
which have an amplitude in excess of a specific value as a function of the amplitude; 

g. means for comparing the at least one curve with predetermined standard cumulative count curves and for 
20 deducting at least one feature from a set of features comprising: material properties, particle concentration, particle 

shapes, particle size and standard deviation thereof and particle size distribution. 

In one embodiment, such equipment is characterized in that 

25 - during operation, the acoustic source generates an acoustic beam which has such a large aperture, and makes the 
time which elapses between two successive measurements so short, that each particle is exposed several times 
while passing through the beam and that, depending on the lateral position of the particles in the measurement vol- 
ume, a different angle-dependent reflection of the acoustic signal is produced; 

identification means are also present for identification of one or more different groups of particles present in the 
30 suspension on the basis of a series of successive angle-dependent reflection signals received over time; 

the means for composing the cumulative curve compose the curve, the maximum value of the amplitudes of a 
series of successive angle-dependent reflection signals, received over time, from a detected particle being taken 
as the amplitude of the electrical signals, which are produced after conversion of the acoustic reflection signals 
from that particle. 

35 

In the last-mentioned equipment, use is advantageously made of the angle-dependent reflection behaviour of 
acoustic signals reflected by particles and it is no longer necessary to meet the condition that only one particle may be 
present in the measurement volume at any one time, as has been explained above. 

In an alternative embodiment the invention relates to equipment for the detection and identification of particles in a 
40 suspension, comprising: 

a. an acoustic source for the generation of acoustic signals; 

b. means for directing the acoustic signals at at least one measurement volume within the flowing suspension, the 
boundaries of the measurement volume in the axial direction with respect to the acoustic source being defined with 

45 the aid of time windows; 

c. means for receiving acoustic reflection signals produced by reflection of the acoustic signals by the particles in 
the at least one measurement volume; 

d. means for converting the acoustic reflection signals into electrical reflection signals; 

e. means for counting numbers of electrical reflection signals which have an amplitude in excess of a predeter- 
so mined value and for converting said count into numbers of particles which are larger than a certain size; 

characterized in that the equipment also comprises: 

f. means for applying an inversion algorithm on the amplitudes of the electrical reflection signals to deduce at least 
one feature from a set of features comprising: material properties, particle concentration, particle shapes, particle 
size and standard deviation thereof and particle size distribution. 

55 

The invention will be explained in more detail below with reference to a few drawings and illustrative embodiments, 
which are intended solely for the purposes of illustration and not as a limitation of the inventive concept. 
In the drawings: 
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Figure 1 shows a diagrammatic measurement set-up for the determination of the particle concentration and particle 
size distribution in a flowing suspension containing equivalent particles; 

Figure 2 shows a measured cumulative count curve for oil droplets in water (133 ppm), for an oil droplet average 
diameter of 40 urn; 

5 Figure 3 shows a diagrammatic measurement set-up for the determination of the particle concentration and particle 

size distribution in a flowing suspension, where several particles may be present in the measurement volume; 
Figure 4 shows simulated cumulative count curves for oil droplets in water, at a concentration of 9.10 9 m' 3 , for a 
standard deviation of the oil droplet diameters of 1 and an average of oil droplet diameter varying from 15 um 
to 20 um; 

10 Figure 5 shows simulated cumulative count curves for oil droplets in water, for a standard deviation of the oil droplet 
diameters of 1 um and for an average of oil droplet diameter of 1 5 um at a concentration varying from 5.1 0 9 m" 3 to 

1l.l0 9 m* 3 . 

Figure 6 shows simulated cumulative count curves for oil droplets in water, at a concentration of 9.10 9 m 3 , for an 
average of oil droplet diameter of 15 ^im and for a standard deviation of the oil droplet diameters varying from 0.5 
15 to 2.5 um; 

Figure 7 shows an inversion result of a inverted simulated count curve. The particle diameter is distributed accord- 
ing to a Gaussian distribution with an average particle diameter u D of 10 urn, a standard deviation a D of 1 jam and 
the particle concentration equals 7.10 9 m" 3 . 

Figure 8 shows an inversion result of a suspension with two particle components: one has a average diameter of 6 
20 jim, the other of 1 0 um. 

Figure 9 shows successive recordings for a single spherical particle in a flowing suspension using the measure- 
ment set-up according to Figure 3; 

Figure 1 0 shows successive recordings of a single elongated particle with a length of 1 mm in a flowing suspension 
using the measurement set-up according to Figure 3; 
25 Figure 1 1 shows a sequence of recordings for two particles 0.5 mm apart. 

First of all a method for the determination of the particle size distribution and particle concentration in a suspension 
containing equivalent particles will be explained. Equivalent particles are understood to be particles of the same mate- 
rial and of the same shape, whilst the size of the particles may vary. 

30 Not only is the number of recordings above the threshold value determined, but a complete histogram of recorded 
amplitudes is made. Said histogram is then converted to give a cumulative count N c with the highest recorded amplitude 
as the start value. The cumulative count N C (A) is defined as the number of recordings with an amplitude greater than or 
equal to A. Said cumulative count curve is simple to calculate from the histogram. The particle concentration and the 
particle size distribution in a suspension containing equivalent particles can be estimated from the histogram or from 

35 the count curve. 

Cumulative count curves for suspensions containing equivalent particles can be simulated on the basis of the 
model given below. First of all it has to be checked whether the assumption that the likelihood of finding two or more 
particles is negligibly small is realistic. 

The Poisson distribution can be used for the probability density function f(m) for finding m particles in the measure- 
40 ment volume: 

'<"»>- < 1 > 

45 

In this equation p is the expectation value for the number of particles in the measurement volume. For suspensions with 
low particle concentrations and for a sufficiently small measurement volume, the likelihood of two particles being 
present in the measurement volume at the same time is negligibly small. In this case it can be deduced from equation 
(1) that the following approximation applies for the likelihood of the presence of no particles (f(0)) and of one particle 
so (f(1)), respectively: 

f{0) = e" p (2) 
55 /(1) = 1-e' p 



The likelihood of finding a particle can also be calculated as follows: 
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'0>- ^measC (3) 



where V meas is the measurement volume and C the number of particles per unit volume. 
Equating f(1) from equations (2) and (3) gives an estimate for p:Sp 



Sp = -m<1 - C) (4) 



Using this equation it is then possible, with the aid of equation (1 ), to make an estimate of the likelihood of finding 
10 two or more particles, so that it is possible to investigate whether the assumption that this likelihood is negligibly small 
is realistic. 

The shape of the cumulative count curve depends not only on the properties of the suspension but also on the 
acoustic pressure distribution in the measurement volume. The probability density function for the measured amplitude 
for a particle of given diameter D depends on the properties of the acoustic field and will be represented by g(A|D). If 
is the suspension consists of a collection of particles of different dimensions, the resultant cumulative count curve will be 
the weighted sum of the individual cumulative count curves which are associated with the particles of a constant diam- 
eter. Said weighting consists in the application of the particle size distribution for the particles in the suspension con- 
cerned. The following equation then applies for the probability density function p(A) for the amplitude: 

20 p{A) = l g{A\D) h(D)dD (5) 



In this equation h(D) is the probability density function for the particle diameter D. 

If said measurement volume is chosen to be around the focus of the acoustic beam 20 (Figure 1), the amplitude of 
25 the incident acoustic pressure can be approximated by a Gaussian function: 

P(r) = e' Kf2 (6) 

30 where 

K = -iy (7) 

32^ 

35 



In this equation r is the lateral distance from the principal axis (see Figure 1), k is the wave number and f N is the f- 
number of the transducer, which is defined as the quotient of the distance from the focus, z f , and th£ transducer diam- 
eter (f N a z t /2a , , where a t is the radius of the transducer). 
40 The likelihood of finding no particles, given in equation (2), indicates which fraction of the total number of measure- 
ments yields an amplitude which is less than or equal to A n (amplitude of the noise). If a particle is positioned that far 
form the central axis of the acoustic beam that the reflection amplitude is below the noise level A n , the particle is said 
to be outside the measurement volume V m9as . As such the amplitude of the noise determines the lateral dimensions of 
the measurement volume; after all, with the aid of equation (6) it follows that: 




50 V meas = n R 2 meas (t 2 -t,)*Cf2 . (9) 

In this equation R mea3 is the radius of the measurement volume in the lateral direction and A max is the maximum ampli- 
tude detected for the suspension (largest particle located on the principle axis of the beam). V meas is the measurement 
55 volume which is delimited in the direction of the principal axis z by the chosen time window [t 1 . t 2 ]. 

For those measurements where there is a particle found in the beam, the measured amplitude will depend on the 
precise location of the particle concerned in the measurement volume and on the dimensions of the particle concerned. 

Assuming a Gaussian acoustic beam profile in the lateral direction (Equation (6)) in the measurement volume and 
negligible amplitude variation in the axial direction within the selected time window, it can be deduced that the condi- 
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tional probability density function g(A|D) for the measured amplitude A originating from an arbitrary particle with diam- 
eter D is given by: 



g(A\D) 



if <A <Aq{D) 



2 1 

A lnA 0 (D) -lnA^ 

0 if A > A 0 {D) A 



(10) 



10 



15 



20 



25 



The following relationship can be derived for the conditional cumulative count G(A > A): 

ln* 0 <£>) - InA 



G(A > A \D) = 



1 if A < An 

0 if A > A 0 {D) 



In these equations Aq is the amplitude which would be detected from the particle concerned if said particle were located 
on the principal axis of the beam. Said Aq depends on the size of the particle. For spherical particles, the diameter of 
which is much smaller than the wavelength, Aq is proportional to the third power of the diameter D (Rayleigh scattering): 

A 0 = yD 3 (12) 



In this equation y is a proportionality constant, which depends both on the transducer properties and material char- 
acteristics of the particle, like density and compressibility. 
30 The probability density function of equation (5) is the so-called forward model of the experiment, from which the 
measured amplitude histogram, given the properties of the suspension (particle diameter distribution h(D)) and the 
properties of the acoustic field, can he calculated. This forward model needs to be inverted to be able to deduce the 
properties of the suspension (h(D)) from the measured amplitude. For a Gaussian pressure field, a method for inversion 
is described below. 

35 The integration in equation (5) is replaced by summation and equations (10) and (12) are applied: 

p<*/>-Zjt— br rt(D / )A0 (13) 

AO "A n ' 

where the summation is over all particle diameters Dj, which might give a reflection amplitude equal to Aj. 

The likelihood for a particle to have a diameter in the interval Dj- AD/2.. Dj+ AD/2 is given by h(Dj)AD, and will be 
45 denoted as q(Dj). So equation (13) becomes: 

p<*/>-Zi— i ^ d j) < 14 > 

so A n ' 

This equation can be written as a matrix product of forward matrix W with probability vector q. 
55 $(A)= W(A,D) 5(D) (15) 

No assumptions were made about the particle diameter distribution, it was only assumed that there exists a maximum 
diameter, which is not a limiting restriction. 
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Equation (15) was deduced for a Gaussian pressure field, however, in general it is possible to deduce a matrix pres- 
entation for the relation between the particle size distribution q(D) and the amplitude distribution p(A). 

There are several ways to find a (stable) pseudoinverse of the forward matrix W which fits the parameters (i.e. the 
elements of vector q) to the data in a least squares sense. A standard approach for this problem is singular value 
5 decomposition (SVD). 

Another, more superior method, uses a stochastic approach, which is described by Franklin (Franklin, J. N., 1970, 
"Well-posed stochastic extension of ill-posed linear problems", Journal of mathematical analysis and applications, vol. 
31 . pp 682-71 6 ( 1 970)). A stochastic approach means that it is assumed that the parameters (q) and noise are samples 
drawn from random processes. The assumption yields in a smoothness constraint for the desired solution. The param- 

10 eters that control the smoothness of the solution depend on the solution itself, therefore an iterative method can be used 
to find the optimal set of parameters. 

The calculated particle size distribution, using one of above-mentioned techniques is not the true size distribution 
of particle in the suspension (q tru e). but the apparent size distribution (q app ). This distribution deviates from the true size 
distribution since smaller particles cannot be detected throughout the complete measurement volume as defined by 

15 equation (9). Applying a correction factor results in the true particle size distribution of particles in the suspension. 
Assuming a Gaussian pressure field, the correction factor is given by: 

,n( ^ ) 

20 Q t rue(D)^^q app (D) (16) 

25 In this equation D max is the diameter of the largest particle in suspension. The factor 4 is applied to normalize the 

area below q true equal to unity. 

The number of measurements with an amplitude larger than A n is indicative for the number of particles C per unit 
volume in the suspension. Assuming a Gaussian pressure field, it can be deduced that the number of particles per unit 
volume is given by: 

30 

c-£t^J«.~.<H <"> 

,w tot v meas J 

35 

In this equation N(...) gives the number of measurements for which the condition given in brackets applies. N tot is 
the total number of measurements; V meas is the measurement volume as defined in equation (9). The factor within 
straight brackets Qisa correction factor which has to be applied because smaller particles cannot be detected within 
the whole measurement volume. 
40 Figure 2 gives an example of a measured cumulative count curve for oil droplets in water (1 33 ppm) , for an oil drop- 
let average diameter of 40 ^m. 

The count curves given in Figures 4 to 6 are the results of simulations based on the combination of equations (5), 
(10) and (12). 

The oil droplets are homogenously distributed throughout the suspension. For the simulations the diameter of the 
45 droplets is normally distributed with an average of jiq and standard deviation ar> The number of droplets per unit vol- 
ume is C. The cumulative count curve for such a suspension depends on these three variables CT D and C. This is 
illustrated in Figures 4, 5 and 6, which show simulated count curves on the basis of the model described above. 

Figure 4 gives simulated count curves for water-oil suspension where n D varies between 15 and 20 fim, where <j d 
is kept constant at 1 |im and C is constant 9.10 9 m" 3 . Figure 5 gives simulated count curves for water-oil suspension 
so where a D is varied from 0.5 to 2.5 urn, the concentration is 9.10 9 m" 3 , and the average drop diameter (\i D ) is 1 5 jim. 

Comparing Figure 4 and 5 show clearly that the effect of a variation in n D differs from the effect of a variation of a D . 
The effect of a variation in C on the cumulative count curve is illustrated in Figure 6. In this Figure, the cumulative count 
curves for a concentration of 5.10 9 m~ 3 , 7.10 9 m" 3 , 9.10 9 m" 3 and 11. 10 9 m" 3 . where m d has been kept constant at 15 
um and a D has been kept constant at 1 jim, are shown. The effect of a higher C is to raise the height of the cumulative 
55 count curve for all A > Ap. Comparing Figure 6 and Figure 4 shows that the effect of a variation in C can be readily dif- 
ferentiated from that of a variation in whilst comparing Figure 6 and Figure 5 shows that the effect of C can also be 
readily differentiated from the effect of or> 

It follows from this example that the C. \x D and o D can be derived from the measured cumulative count curve for a 
suspension in which an unknown quantity of oil droplets of unknown ^ D and a D are distributed. For this purpose meas- 
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ured cumulative count curves can be compared with predetermined standard cumulative count curves for known parti- 
cles. Said standard cumulative count curves can, for example, be determined via simulation, as has been explained 
above with reference to Figures 4, 5 and 6. In this case, for example, a Bayesian inversion can be used to derive the 
particle size distribution from the cumulative count curves. As an alternative, standard cumulative count curves of this 

5 type can be determined experimentally. Standard count curves of this type can be stored in the memory 27 of computer 
26, the computer 26 being equipped with suitable software for comparing measured count curves with the standard 
count curves. Standard methods, such as the least squares fit. can be used for this purpose. The method can, of 
course, also be used for suspensions containing particles which have a particle size distribution other than normal par- 
ticle size distribution with characteristic values other than n D and a D . 

w Another method to derive C and the particle size distribution (e.g. defined by n D and a D ) from a measured cumu- 
lative count curve is to apply an inversion algorithm, such like the one as given in equation (1 5). The capabilities of such 
an inversion algorithm are illustrated in Figures 7 and 8. 

Figure 7 gives the inversion result using the described inversion technique given by equations (14), (15), (16). The 
line (curve 16) shows the exact size distribution, the circular markers (curve 17) are the result obtain using singular 

15 value decomposition and the cross markers (curve 18) are the inversion result using the stochastic approach. The Fig- 
ure shows that more accurate results are obtained using the stochastic approach. To illustrate that any kind of particle 
size distribution can be determined using this technique, an inversion result of a non-gaussian size distribution is shown 
in Figure 8. The line (curve 19) gives the exact distribution and the cross markers (curve 20) give the inversion result 
using the stochastic approach. 

20 A method and set-up for the characterisation of various types of particles in a suspension will now be described. 

In the above it has been assumed that there is only one type of particle in the suspension, only the size of the par- 
ticles being allowed to vary. A further assumption was that the likelihood of more than one particle being present in the 
measurement volume at the same time is negligibly small. On these grounds, a small measurement volume in the focus 
region of the acoustic beam 20 was preferably chosen. 

25 However, there are sometimes advantages in, in contrast, selecting a large measurement volume beyond the 
focus. This is illustrated in Figure 3. in which the same reference numerals as in Figure 1 refer to the same elements. 
The measurement volume is delimited in the axial direction of the principal axis z by z 3 and z 4 associated with time win- 
dow [t 3 . U]. In this case the angle of incidence within the measurement volume varies as a function of the lateral position 
r with respect to the principal axis z. If several reflection measurements are carried out in succession on one particle 22 

30 while the particle is passing through the measurement volume (see Figure 3), the angle-dependent reflection behaviour 
of the particle can be determined from the change in the reflection signal as a function of the lateral position. Said 
angle-dependent behaviour is highly dependent on the shape of the particle and therefore the particle can be charac- 
terized on the basis of this behaviour. 

Figure 9 shows a succession of simulated recordings for a spherical particle which passes through the beam at a 

35 fixed distance beyond the focus. Simulated recordings for an elongated particle are given in Figure 10. These Figures 
clearly illustrate the effect of the shape of the particle on the reflection behaviour as a function of the angle. 

For this method, the larger the aperture of the beam from the transducer 23. the more clearly visible is the differ- 
ence in angle-dependent behaviour between different types of particles. As a result, the measurement volume therefore 
also increases. However, the requirement that no more than one particle at a time may be present in the lateral direction 

40. in the measurement volume is dispensed with with this method because different particles located laterally alongside 
one another can be individually recognised in the succession of recordings. This is illustrated in Figure 1 1 , where two 
spherical particles pass through the beam just behind one another. Although the lateral distance (15 mm) between 
these particles is much smaller than the beam diameter, the reflections originating from the first particle can be sepa- 
rated from the reflections from the second particle in the succession of recordings. The most important characteristic 

45 of the last-mentioned embodiment is that the beam has a large aperture and that the time which elapses between two 
successive measurements is so short that each particle is exposed several times during its presence in the beam. As 
has been explained above, this can be achieved, for example, with the aid of a focused transducer. Theoretically, the 
use of a point source or a relatively small source (of the order of magnitude of the wavelength used) is also possible. It 
is only at high frequencies that this cannot be achieved in practice. Optionally, use can be made of a Radon transfor- 

50 mation or a 2D Fourier transformation to the k x -k domain in order to quantify the angle-dependent reflection behaviour 
of each recorded particle. 

Following the characterisation of a particle on the basis of its angle-dependent behaviour, the cumulative count 
curve can be determined for each type of particle. For this purpose the maximum amplitude for the series of recordings 
originating from the particle concerned (this is the amplitude when the particle is located on the principal axis z Of the 
55 beam 20) is taken as the amplitude per detected particle. This maximum amplitude can still vary for one type of particle 
because of variations in particle size and variations in the minimum lateral distance from the principal axis z at which 
the particle passes through the beam. A cumulative count curve similar to that illustrated above with reference to Fig- 
ures 4, 5 and 6 will be found for each type of particle. 

A possible application of the method described above is the measurement of thrombus particles in a blood stream. 
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Patients who have a heart valve prothesis have an increased risk of thrombosis because the artificial valve promotes 
the production of thrombus particles. Therefore, thes patients have an increased risk of acute vascular occlusion as a 
result of too large a thrombus particle. This can lead, for example, to cerebral infarction. 

Anticoagulants are administered to the patient to counteract thrombus formation. The risk of too high a dosage of 

5 these anticoagulants is the occurrence of haemorrhaging (for example cerebral haemorrhaging). Currently the dosage 
is determined on the basis of the coagulation measured on a blood sample. However, this is an unreliable method 
because the coagulation is dependent on many more factors than solely the concentration of thrombi. Moreover, a 
measurement of this type provides only a snapshot. 

It is therefore desirable to have available a method with which the concentration of thrombi in the blood can be 

10 measur d reliably. The measurement method described here is suitable for this purpose. 

By using a transmission frequency of the sound waves which is relatively high for medical applications, for example 
in the range of 1 0-40 MHz, preferably 20-30 MHz, it is possible to detect small particles present in the blood. The depth 
of penetration is still found to be sufficiently great at these frequencies. By making use of the high efficiency of, for 
example, a composite transducer and of sophisticated transmission and reception electronics, it is possible to achieve 

is an adequate depth of penetration, so that an echo-acoustic recording of the blood can be made using a non-invasive 
technique. However, blood also contains other particles which will be detected by the ultrasonic reflection method, for 
example red blood cells. 

It is known that red blood cells form elongated aggregates during a certain period of the heart cycle (end of diasto- 
lic), which aggregates align with the flow. See, for example, M.G.M. de Kroon: "Acoustic backscatter in arteries - Meas- 
20 urements and modelling of arterial wall and blood", thesis 1993. ISBN 90-9006182, Section III. It is possible to 
distinguish these long structures from thrombi, which are of a jagged shape (see: S. Chien: "Clinical Haemorheology", 
Martinus Nijhoff). on the basis of the angle-dependent behaviour. 

After making this distinction, the concentration and size distribution of thrombi can be estimated on the basis of the 
cumulative count curves, as has been explained above. 

25 

Claims 

1 . Method for the detection and identification of particles in a suspension, comprising the following steps: 

30 a. generation of acoustic signals having the form of a beam using an acoustic source; 

b. directing the acoustic signals at at least one measurement volume within the suspension, the boundaries of 
the measurement volume in the axial direction with respect to the acoustic source being defined with the aid of 
time windows; 

c. reception of acoustic reflection signals produced by reflection of the acoustic signals by the particles in the 
35 at least one measurement volume; 

d. conversion of the acoustic reflection signals into electrical reflection signals; 

e. counting numbers of electrical reflection signals which have an amplitude in excess of a predetermined 
value and conversion thereof into numbers of particles which are larger than a certain size; 
characterized in that the method also comprises the following steps: 

40 f. composing at least one curve on the basis of a cumulative count of the number of reflection signals which 

have an amplitude in excess of a specific value as a function of the amplitude; 

g. comparison of the at least one curve with predetermined standard cumulative count curves and deduction 
of at least one feature from a set of features comprising: material properties, particle concentration, particle 
shapes, particle size and standard deviation thereof and particle size distribution. 

45 

2. Method according to Claim 1 , characterized in that the standard cumulative count curves are derived as follows: 

definition of the lateral boundaries of the measurement volume as being the boundaries beyond which no 
amplitudes of the reflection signals smaller than a predetermined minimum value can be recorded; 
so - derivation of a probability density function for the amplitude of the electrical ref lection signals within the meas- 

urement volume as a function of a given particle shape, and particle size distribution thereof, and 
derivation of the standard count curves from the probability density function for the amplitude. 

3. Method according to Claim 2, characterized in that the following equation is used for the probability density function 
55 for the amplitude: 

P(A) = I g(A\D) h(D)dD 

where: 
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4. 



g(A|D) = probability density function for the measured amplitude A originating from an arbitrary particle with 
diameter D; 

h(D) = probability density function for the particle diameter, 

the equation being integrated between two predetermined limits for the particle diameter, for example 0 and «>. 

Method according to Claim 3, characterized in that a Gaussian amplitude profile of the acoustic signals in the lateral 
direction of the measurement volume and a substantially negligible amplitude variation in the acoustic signals in the 
axial direction of the measurement volume are assumed for determination of g(A|D), resulting in the relation: 



g(A\D) 



if A^ <A <A 0 (D) 



1 1 

~A lnA 0 (X» - ln^ 

0 if A > Aq(D) A A<A^ 



15 



where Aq is the amplitude which would be detected from the particle concerned if said particle were located on the 
principal axis of the beam. 

20 

5. Method according to Claim 1, characterized in that standard cumulative count curves are determined experimen- 
tally with the aid of suspensions containing particles of a pre-known shape, size and standard deviation of the size. 

6. Method for the detection and identification of particles in a suspension, comprising the following steps: 

25 

a. generation of acoustic signals using an acoustic source: 

b. directing the acoustic signals at at least one measurement volume within the suspension, the boundaries of 
the measurement volume in the axial direction with respect to the acoustic source being defined with the aid of 
time windows; 

30 c. reception of acoustic reflection signals produced by reflection of the acoustic signals by the particles in the 

at least one measurement volume: 

d. conversion of the acoustic reflection signals into electrical reflection signals; 

e. counting numbers of electrical reflection signals which have an amplitude in excess of a predetermined 
value and conversion thereof into numbers of particles which are larger than a certain size; 

35 characterized in that the method also comprises the step of applying an inversion algorithm on the amplitudes 

of the electrical reflection signals to deduce at least one feature from a set of features comprising: material 
properties, particle concentration, particle shapes, particle size and standard deviation thereof and particle 
size distribution. 

40 7. Method according to Claim 6, characterized in that the following equation is used for the inversion algorithm: 

p(A)= W{A,D)q(D) 

45 where W is a forward matrix of size j x i, i and j being predetermined integer numbers, which relates likelihood q(Dj) 

for a particle to have a diameter in a predetermined interval Dj-AD/2..Dj+AD/2 to probability p(Aj) of measuring an 
amplitude A,. 

8. Method according to Claim 7, where the forward matrix W is inverted applying either a Singular Value Decomposi- 
50 tion or a stochastic approach. 

9. Method according to claim 6, characterized in that a Gaussian amplitude profile of the acoustic signals in the lateral 
direction of the measurement volume and a substantially amplitude variation in the acoustic signals in the axial 
direction of the measurement volume are assumed for determination of the forward matrix W. resulting in the fol- 

55 lowing expression for p(A,): 
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1 0. Method according to claim 6, where a calculated particle size distribution resulting from the inversion algorithm 
is corrected for the fact that smaller particles cannot be detected throughout the whole measurement volume, 
resulting in a true particle size distribution q true . 

70 

11. Method according to claim 10, characterized in that a Gaussian amplitude profile of the acoustic signals in the lat- 
eral direction of the measurement volume and a substantially negligible amplitude variation in the acoustic signals 
in the axial direction of the measurement volume are assumed for determination of the said correction to be 
applied, resulting in the following relation: 



75 



Y°max 



20 y D 3 



where D max is the diameter of the largest particle in suspension and £ is a factor applied to normalize the area 
below q true equal to unity. 

25 

1 2. Method according to claim 1 where the particle concentration is calculated from the number of measurements with 
amplitude above the noise level. 

1 3. Method according to claim 12, where the calculated particle concentration is corrected for the fact that smaller par- 
30 tides cannot be detected through the whole measurement volume, resulting in a true particle concentration. 

14. Method according to claim 13, characterized in that a Gaussian amplitude profile of the acoustic signals in the lat- 
eral direction of the measurement volume and a substantially negligible amplitude variation in the acoustic signals 
in the axial direction of the measurement volume are assumed for determination of the said correction to be 

35 applied, resulting in the following relation for the true particle concentration C: 

' v tot u meas |_S ^ J 

40 



where N(...) gives the number of measurements for which the condition given in brackets applies, N tot is the total 
number of measurements, V meas is the measurement volume, £ is a normalizing factor and q true is a true particle 
size distribution. 

45 

15. Method according to one of the preceding claims, characterized in that the measurement volume is chosen to be 
around the focus of the acoustic beam and to be so small that the likelihood of the presence of more than one par- 
ticle in the measurement volume during a measurement is negligibly small. 

so 16. Method according to one of Claims 1 to 5, characterized in that: 

in step a, the acoustic beam generated has such a large aperture, and the time which elapses between two 
successive measurements is so short, that each particle is exposed several times while passing through the 
beam and that, depending on the lateral position of a particle in the measurement volume, a varying arigle- 
55 dependent reflection of the acoustic signal is produced; 

prior to step f. one or more different types of particles present in the suspension are identified on the basis of 
a series of angle-dependent reflection signals received successively over time; 

when composing the curve in step f. the maximum value of the amplitudes of a series of successive angle- 
dependent reflection signals, received over time, from a detected particle from the one or more groups is taken 
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as the amplitude of the electrical signals, which are produced after conversion of the acoustic reflection signals 
from that particle. 

17. Method according to Claim 16, characterized in that in step a. a focused acoustic beam with a predetermined focus 
s is generated and in step b the focus remains outside the measurement volume. 

18. Method according to Claim 16, characterized in that an acoustic source having dimensions smaller than the wave- 
lengths of the generated acoustic signals is used. 

w 1 9. Method according to one of the preceding claims, characterized in that the acoustic signal comprises sound signals 
having a frequency of 10-40 MHz, preferably 20-30 MHz. 

20. Equipment for the detection and identification of particles in a suspension, comprising: 

75 a. an acoustic source for the generation of acoustic signals; 

b. means for directing the acoustic signals at at least one measurement volume within the flowing suspension, 
the boundaries of the measurement volume in the axial direction with respect to the acoustic source being 
defined with the aid of time windows; 

c. means for receiving acoustic reflection signals produced by reflection of the acoustic signals by the particles 
20 in the at least one measurement volume; 

d. means for converting the acoustic reflection signals into electrical reflection signals; 

e. means for counting numbers of electrical reflection signals which have an amplitude in excess of a predeter- 
mined value and for converting said count into numbers of particles which are larger than a certain size; 
characterized in that the equipment also comprises: 

25 f. means for composing at least one curve on the basis of a cumulative count of the number of reflection signals 

which have an amplitude in excess of a specific value as a function of the amplitude; 

g. means for comparing the at least one curve with predetermined standard cumulative count curves and for 
deducing particle concentration and/or particle characteristics, such as particle size distribution, particle shape 
and standard deviation thereof, from said comparison. 



30 



21. Equipment for the detection and identification of particles in a suspension, comprising: 



a. an acoustic source for the generation of acoustic signals; 

b. means for directing the acoustic signals at at least one measurement volume within the flowing suspension, 
35 the boundaries of the measurement volume in the axial direction with respect to the acoustic source being 

defined with the aid of time windows; 

c. means for receiving acoustic reflection signals produced by reflection of the acoustic signals by the particles 
in the at least one measurement volume; 

d. means for converting the acoustic reflection signals into electrical reflection signals; 

40 e. means for counting numbers of electrical reflection signals which have an amplitude in excess of a predeter- 

mined value and for converting said count into numbers of particles which are larger than a certain size; 
characterized in that the equipment also comprises: 

f. means for applying an inversion algorithm on the amplitudes of the electrical reflection signals to deduce at 
least one feature from a set of features comprising: material properties, particle concentration, particle shapes, 
45 particle size and standard deviation thereof and particle size distribution. 

22. Equipment according to Claim 20, characterized in that 

during operation, the acoustic source generates an acoustic beam which has such a large aperture, and 
so makes the time which elapses between two successive measurements so short, that each particle is exposed 

several times while passing through the beam and that, depending on the lateral position of the particles in the 
measurement volume, a different angle-dependent reflection of the acoustic signal is produced; 
identification means are also present for identification of one or more different groups of particles present in the 
suspension on the basis of a series of successive angle-dependent reflection signals received over time; 
55 - the means for composing the cumulative curve compose the curve, the maximum value of the amplitudes of a 
series of successive angle-dependent reflection signals, received over time, from a detected particle being 
taken as the amplitude of the electrical signals, which are produced after conversion of the acoustic reflection 
signals from that particle. 
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fig-4 




amplitude [V] 



C = 9.10* m 3 ; 
O d = 1 nm; 

1: Md = 15 
2: ji D = 16 p.m; 
3: n D = 17 ^m; 
4: |i D = 18 jim; 
5 \i D = 19 nm; 
6: \i D = 20 jim, 
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C = 9.10 9 m' 3 ; 
H D = 15 urn; 
7: C D = 0.5 urn; 
8: a D = 1.0 urn, 
9: a D = 1.5 nm, 
10: a D = 2.0 iim\ 
11: a D = 2.5 prn; 
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fig-b 




amplitude [V] 

H D = 15 Mm; 
a D = 1 |am, 
12: C=5.10 9 m" 3 
13: C=7.10 9 m- 3 
14: C=9.10 9 m 3 
15: Oil. 10* m*' 
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16: real distribution ( ) 

17: inversion result with Singular Value Decomposition (o) 
18: Inversion result with stochastic approach (+) 
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